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ION  BEAM  PROPAGATION  IN  A  FILAMENTED  CHANNEL 


I.  INTRODUCTION 

The  production  of  intense  focused  ion  beams ^  has  led  to  the 
consideration  of  using  z-discharge  plasma  channels  to  transport  them 
several  meters  to  inertial  fusion  targets?.  Analysis  has  shown  that 
although  electrostatic  beam-plasma  streaming  modes  are  stable-*,  electro¬ 
magnetic  microinstabilities  will  occur  during  transport  in  such 
channels^.  The  fastest  growing  mode  produces  current  bunching  of  the 
electrons  in  the  channel.  In  this  report  ion  beam  propagation  in  a 
filamented  channel  is  investigated  in  order  to  determine  the  effects  on 
radial  beam  containment  in  channels  and  on  the  radial  beam  density 
profile. 

In  Sec.  II,  radial  current  bunching  in  the  channel  is  considered 
and  in  Sec.  Ill,  azimuthal  current  bunching  is  considered.  Finally 
the  results  of  this  work  are  summarized  in  Sec.  IV. 


II.  RADIAL  CURRENT  BUNCHING  IN  THE  CHANNEL 

If  k  =  k  er>  then  radial  current  bunching  occurs  in  the  channel 
and  the  net  current  density  can  be  modeled  by 

i  =[l(r)  +  j1  H(r-rb)  JQ(kr)j  §z  ,  (1) 


where 


H(r-rb) 


0 


S  r  s;  r, 
D 

r  >rb 


i 
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r,  is  the  beam  radius  and  J  is  the  zero  order  Bessel  function 
b  o 

which  satisfies  the  wave  equation  in  cylindrical  coordinates.  For 
Note:  Manuscript  submitted  January  2,  1980. 


a  channel  of  radius  r 

c 

r  r 

1  =  2  it  f C  j(r)  rdr  =  2  n  fC  j  (r)  rdr 

*o  *0 

is  the  channel  current  which  confines  the  ion  beam.  Here  k  =  a.0/r, 
was  used,  where  a^  is  the  zero  of  first  order  Bessel  function 
J^.  The  second  term  on  the  right  hand  side  of  Eq.  (1)  is  then  the 
residual  current  density  arising  from  the  bunched  electron  current 
superimposed  on  the  unbunched  ion  beam  current.  Since  j (r)  is 
established  in  an  initially  cold  plasma  and  is  driven  by  a  capacitor 
bank  on  a  time  scale  much  longer  than  the  beam  pulse  duration, when 
the  beam  is  injected  into  the  now  highly  conducting  plasma  the  total 
net  current  resists  change.  Thus,  in  the  absence  of  bunching 
j  =  je  +  jb  =  j<r>  6a,  where  j  «  jfe. 

The  magnetic  field  which  determines  the  ion  motion  in  a  filamented 
channel  is  then  found  from 


1  d 
r 


or  (r  V  =  f  [j(r)  +  h  Jo(kr)]  ’rSrt 


(2) 


—  S”1 

Assuming  j (r)  =  j  (r/r^)  inside  the  channel  (s  a  1) , 


Be  = 


4n  jQ  r- 


(s+1)  c  r 


s-1 


4n  j. 


kc 


(kr)  ,  r  £  r^  , 


(3) 


where  r^  >  r^.  The  parameter  s  is  a  magnetic  field  shaping  factor  which 
is  equal  to  one  at  the  front  of  the  beam  and  increases  as  one  moves 
toward  the  tail  of  the  beam.  The  beam  heated  plasma  in  the  channel  at 
the  tail  of  the  beam  expands  radially  carrying  the  channel  current  with 
it.  This  decreases  the  magnetic  field  (increases  s)  near  the  center  of 
the  beam  and  produces  a  sharply  rising  magnetic  field  at  the  edge  of  the 
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beam.^ 

The  total  energy  H,  the  axial  cononlcal  momentum  P  ,  and  the  canon- 

z 

ical  angular  momentum  Pq  are  all  constants  of  the  motion  for  the  ions. 
From  Az  *  "  f  B^dr  aru*  ®q*  (3),  P^  is  found  to  be 


2uicbrc 

(8+1) 


/  \s+l  2uu  ,  (s+l)  j  T  "I 

(i)  -  ■ 

'  f  co  J 


2 

where  uu  ,  =  2ner  1  /m.  c  (s+l)  .  When  j./j  >  0,  it  is  clear  that  the 
cb  co  i  1  o 

s+l 

beam  ion  orbits  are  radially  confined  since  1  -  J  (kr)  £  0  and  r  is 

o 

monotonically  increasing.  When  j^/j  <  0,  the  ion  orbits  are  also  radi¬ 


ally  confined  if 


is  sufficiently  small  since 


[l  -  Jo(kr)J 


1.4  and  r  is  a  monotonically  increasing  function  of 

radius . 

Now  consider  this  result  in  more  detail.  Since  Pq  =  0^,  the  minimum 
axial  canonical  momentum  is  given  by 


/  cos  a  - 
o  o 


2Vc  /i\S+1 

(s+1)  \rc)  J  • 


when  j^  =  0.  Here  a.Q  is  the  maximum  ion  injection  angle,  Vq  =  (2H/m^) 
and  rg  is  the  beam  focal  spot  size  at  the  point  of  injection  into  the 
channel.  The  beam  envelop  radius  is  then 

r, ...  /  p° \  i 8+r 


Is+1L  (v 

2uJcbrc  \°  mi  / 


where  r  *  r°  when  z  =  V  for  an  ion  injected  with  P  *  P°.  When  s  *  1 
b  o  z  z 

the  usual  result  is  recovered^ 


-  ('-■$-)] 

r  V  m.V  / 

3  C  \  i  0/ 


3 


When  j  +  0,  the  minimum  axial  canonical  momentum  is  given  by 


=  P  - 
z  z 


2mjj  ih)cb(s+1> 

j  k2r 
o  c 


[l  -  Vkrs>] . 


(8) 


where  this  result  only  strictly  applies  for  |j^/jQ|  <  krc(rg/rc)S/(s+l) , 
In  this  case  the  beam  envelope  radius  is  found  from 


(s+1) 


/r  \  S+1  21,  (s+1)  r  TV  /  P1  \ 

(ir)  +  “Try  f1-  -7) 

\c/  jkrL  J  ID  ,  r  l  m.V  / 

'  '  Joc  cb  c  \  i  0/ 


(9) 


For 


Vjo 


<<1  this  reduces  to 


H  £?©"'[- 


J  (kr  )  -  J  (kr  ) 
ob  os 


■]) 


(10) 


so  that  for  <  0  the  beam  radius  expands,  whereas  for  j  >  0  the  beam 
radius  decreases.  For  larger  J  j^/jQ  |  Eq.  (9)  must  be  solved  numerically. 

When  j ^  <  0,  beam  ions  can  be  magnetically  trapped  and  prevented 
from  reaching  the  axis  if  J  j j q  J  is  sufficiently  large.  Thus  the  density 
profile  could  be  depressed  on  axis.  For  a  given  ion  injected  into  the 
channel  at  a  radius  rQ  and  at  an  injection  angle  a  trapping  occurs  if 


1  — 

> 

.2  2 
k  r 

c 

V  (1-cos  a)  (r  /r  )s+1 

O  j  0  c 

1  jo 

1-J  (kr  ) 

0  0 

2u),r  (s+1)  (s+1)2 

CD  C 

(11) 


When  s  =  1  and  a  <<  1,  this  reduces  to 


|£i 

.2  2// 
k  r  /4 
v.  c' 

2  2 

V  a  /r  \ 

°  1  (  0  \ 

k 

[1-Jo<kro>] 

.2u,cbrc  \rc/  . 

(12) 


so  that  trapping  begins  when  the  bunched  current  density  exceeds  the  mag¬ 
nitude  of  the  channel  current  density  which  is  required  to  confine  the 
beam.  Furthermore,  Eq.  (9)  shows  that  the  equilibrium  beam  radius  is 
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also  significantly  modified  if  |  j^/j  |>  1. 

Thus  when  radial  current  bunching  occurs  in  the  channel,  good  beam 
propagation  is  still  expected  unless  j  j^/jo  j  exceeds  unity.  If  |j^/jQ| 

>  1,  the  beam  will  be  expelled  from  the  center  of  the  channel  when  j^ 

<  0  or  the  beam  will  pinch  when  j^  >  0.  If  r^  becomes  larger  than  r£ 
the  analysis  breaks  down  and  a  more  sophisticated  analysis  is  required, 
however,  it  has  been  established  here  that  it  is  desirable  to  maintain 

|  j  i/ |  <  **  For  ^i**  Jn  exp  ^Tb^*  8°0<^  keam  Propagation  requires 

j  <  j  exp  (v  t.  ).  Here  v  is  the  growth  rate  of  the  instability,  T, 

nob  •  b 

is  the  beam  pulse  duration  and  j  is  the  initial  level  of  noise  in  the 
electron  current.  Since  the  predominate  source  of  j  is  the  nonuniform¬ 
ities  in  the  radial  profile  of  the  ion  beam  current  density  (remember 
j  -«  j  +  j  ),  jitter  in  the  ion  beam  current  density  must  be  kept  below 

“  <  T1  exp  (-  YTh).  (13) 

Jb  Jb 

4 

Typically  Y ~  1.0  -  2.0  and  j  /j^  ~  0.1  for  proposed  fusion  systems. 

Thus  the  jitter  in  the  beam  current  density  must  be  kept  below  1  -  4%. 

4 

Since  the  filamentation  instability  is  nonconvective  ,  these  results  only 

0 

apply  at  the  tail  of  the  beam.  At  the  front  of  the  beam  the  limit  on 
^bn^b  statec*  in  Eq.  (13)  is  less  severe  and  can  be  found  by  replacing 
y  by  Y^x/xjj  where  x  is  the  distance  from  the  front  of  the  beam  for  a 
beam  of  length  x^. 

HI .  AZIMUTHAL  CURRENT  BUNCHING  IN  THE  CHANNEL 

A 

If  ■  k«g  then  azimuthal  current  bunching  occurs  in  the  channel 
and  the  net  current  density  can  be  modeled  by 
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2  =  (j (r)  +j1H(r-rb)sin  ma)ez  , 


(14) 


where  m  =  3,4,5,....  will  be  considered.  Special  treatment  is  required 
for  m  =  1,2  which  is  not  considered  here  since  m~  3-7  in  typical  cases. ‘ 
Again 


I  = 
o 


/ 


j (r ,0 )  rdr 


■ 2’/ J 


(r)r  dr 


is  the  channel  current  which  confines  the  ion  beam  and  j  sin  m6  is  the 
residual  current  density  arising  from  the  bunched  electron  current  super¬ 
imposed  on  the  unbunched  ion  beam  current.  The  magnetic  field  which 

determines  the  ion  motion  in  a  filamented  channel  is  then 
s 


B, 


4ni  r 
Jo 


8rtj1r 


0  ,  ...  s— 1 

(s+l)cr 

c 


B  * 


4rrj1i 


r  c(m2-4) 


c(m  -4) 

cos  m8  , 


sin  m0  ,  r  £  r,  , 

D 


r  £  r 


(15) 


(16) 


-S-l 


where  again  j (r)  =  j  (r/r  )  was  assumed. 

o  c 

In  this  case  H  and  P  are  constants  of  the  motion  but  P„  is  not. 

z  0 

From  Eq.  (15)  P  is  found  to  be 
z 


•  2,%rc  (l  \ 

(5+1) 2  Vc/ 


s+1  „  2 . 

2uu  cbr  j-^sin  m0 

+  s  ,  r  s  r  . 

rc(m  -4)jo 


(17) 


Solving  Eq.  (17)  for  z  shows  that  as  long  as 

I  j  i  1 


2  . 
m  -4 


(s+1)' 


(18) 


where  r  ^  r,  ,  all  ion  orbits  are  confined  within  the  channel.  For 
c  b 


larger 


|VM 


ions  can  escape  from  the  channel  and  the  beam  density  is 
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gradually  depleted. 


Even  for  small  | j  J Q  |  beam  density  on  axis  decreases  in  time 
since  |  Pg  |  increases  in  time  for  most  ions.  In  fact  Pg  is  constant  only 
for  those  few  ions  which  are  injected  at  0q  =  nn/m  (n  *  1,2,3, ...  ,2m) . 

ror|VJo|«  1 


(IT  8  '  )'  =  -  R(R"-  R8 


t2  /2mj.jCos  m8\ 

\  V"2-4>  / 


(19) 


2  2 

where  R  =  Kr,  T  =  0  t ,  K  «*  m.  UJ  ,/P  r  ,  U  =  P  U)  .  /m.r  and  a  prime 
*  icbzc  zcbic 

signifies  differentiation  with  respect  to  T.  Thus,  to  zero  order  in  the 
small  parameter  | j  J 0  j 

(R20')  »  K2Pq/Q  =  const.  (20) 

o  U 

If  PQ  is  initially  zero  and  R  «  A  cos  (t  +  $)  (implying  small  injection 
d  o 

angles,**  i.e.  a  «  1)  the  time  averaged  increase  in  |Pg  |  can  be  expressed 


as 


<  |PQ  I : 

r  V  a 
o  o 


=  t/T, 


(21) 


where  r  V  a  r(z=o)r(z=o)  and 
o  o 


roa(m  -A)  |JQ/j1| 

mV  cos  m0 
o  o 


(•••#)" 


(22) 


Thus  in  a  time  T,  |Pgj  increases  from  zero  to  r^V^a..  For  the  average  ion 
<  Icos  m0  |>  =  2/n,  r  ~  r,  /2,  a  »  a  /2  and 

^  I  O  9  O  D  O 


TTr.a  (m  -4) 

D  O 


2mV 


(23) 
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Typically  r  ~  0.4  cm,  a  ~  0.2  rad,  m~  5,  V  ~  3.1  *  1CT  cm/sec  and 
do  o 


8  _  i 

uu  ,  «*  4  x  10  sec  so  that 
cb 


2.3  x  10 


-9 


Vjl 


sec. 


(24) 


Thus  the  beam  hollows  out  on  a  time  scale  of  t  even  for  small 
|j^/jo|  ,  although  the  beam  is  still  confined  within  the  channel.  In 
order  to  prevent  this  hollowing  out  from  occurring  during  beam  propagation 
one  needs  T  <  T  or 


nr  a  (m2-4)j  /  „  r2ui  ' \  1 

Jn  <  - 2mL -  K  + 


exp  (~YTb) 


(25) 


where  =  L/V^  is  the  beam  transit  time  in  a  channel  of  length  L  and 
again  j  is  the  initial  noise  level  in  the  electron  current.  Thus  as 
argued  at  the  end  of  Sec.  II  this  implies  that  the  jitter  in  the  ion  beam 


current  density  must  be  kept  below 


fbn  nVo(m  ~4) 
jb  2mL 


^  (  2+  VcbV1 

5T  V" 


exp  (- yt  ) 


,-4 


(26) 


Typically  this  is  found  to  be  j b  ~  1  ~  5  x  10  ’  for  proposed  fusion 
systems.  Since  this  low  level  of  noise  is  probably  not  achievable,  some 
hollowing  out  of  the  beam  is  likely  if  current  bunching  occurs,  however 
the  beam  will  still  be  confined  in  the  channel  for  small  | j  J 0 1  •  This 
effect  will  be  most  pronounced  at  the  tail  of  the  beam  where  largest 

growth  of  the  nonconvective  filamentation  instability  occurs. 

II  2  2 

j^/j  I  ^  (m  -4)/(s+l)  the  beam  density  is  gradually  depleted 
as  increasingly  more  ions  become  unconfined.  This  depletion  will  occur 
in  addition  to  the  hollowing  out  of  the  beam  and  also  occurs  predominately 
at  the  tail  of  beam.  In  order  to  prevent  this  beam  density  depletion 
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one  needs 


Jbn 


(m  -4) 
(s+1)2 


—  exp  (-  YT  b  ) 
Jb 


For  m  ~  5,  s  ~  4,  j  /  j ,  ~  0.1  and  yt , 
ob  b 


(27) 


1  -  2  at  the  tail  of  the  beam. 


•^bn^b  ^  1  "  4%  is  required.  At  the  front  of  the  beam  the  condition  is 
much  less  severe  with  yi^  replaced  by  yi^x/x^  where  x  is  the  distance 
from  the  front  of  the  beam  for  a  beam  of  length  x^.  Eq.  (26)  can  be 
similarly  modified  for  the  front  portion  of  the  beam. 

IV.  CONCLUSIONS 

From  the  analysis  presented  here  it  can  be  concluded  that  good  beam 
transport  in  a  filamented  E-discharge  channel  is  possible  as  long  as  the 
current  bunching  remains  below  certain  levels.  For  radial  current  bunch¬ 
ing  it  was  found  that  should  not  exceed  unity,  which  implies  that 

exp  (-yT^).  If  exceeds  unity,  the  beam  will  be  expelled 

from  the  center  of  the  channel  when  <  0  or  the  beam  will  pinch  when 

h  >  °- 


For  aximuthal  current  bunching  the  beam  hollows  out  on  a  time  scale 

T  defined  in  Eq.  (23).  This  hollowing  out  occurs  for  all  values  of 

11,/j  |  but  will  not  reach  a  significant  level  if  T  <  f  (see  Eq.  (25). 

1  °  *2  C 

For  ^  — — - -  the  beam  density  is  also  gradually  depleted  as  in- 

»  1  °l  (s+1) 

creasingly  more  ions  become  unconfined. 

These  results  set  limits  on  the  level  of  which  can  be  tolerated 
without  seriously  affecting  beam  transport  due  to  current  bunching  effects. 
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